Topological crossovers near a quantum critical point 
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We study the temperature evolution of the single-particle spectrum e(p) and quasiparticle mo- 
mentum distribution n(p) of homogeneous strongly correlated Fermi systems beyond a point where 
the necessary condition for stability of the Landau state is violated, and the Fermi surface becomes 
multi-connected by virtue of a topological crossover. Attention is focused on the different non- 
Fermi-liquid temperature regimes experienced by a phase exhibiting a single additional hole pocket 
compared with the conventional Landau state. A critical experiment is proposed to elucidate the 
origin of NFL behavior in dense films of liquid 3 He. 
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The study of non-Fermi-liquid (NFL) behavior of 
strongly correlated Fermi systems in the regime of a 
quantum critical point (QCP) is currently one of the 
most active and challenging areas of condensed mat- 
ter physicsi 2 - As a rule, such behavior is attributed to 
second-order phase transitions, and the QCP is identi- 
fied with the end point of a corresponding line of tran- 
sition temperatures, denoted by Tm{H) in the proto- 
type in which an external magnetic field H is the con- 
trol parameter. In this case, NFL behavior is triggered 
by critical antiferro- or ferromagnetic fluctuations, which 
lead to violation of respective Pomeranchuk stability con- 
ditions (PSC). Ensuing NFL phenomena are presum- 
ably explained either within the Hertz-Millis theoryM 
or, in heavy-fermion metals, within a Kondo breakdown 
modeljia 2 -^ 

However, the widely promulgated fluctuation scenario 
is inconsistent with experimental data on a number of 
strongly correlated Fermi systems exhibiting NFL behav- 
ior: 

(i) In dense 3 He films where the emergent NFL be- 
havior has been documented, experiments - — has 
not identified any related second-order phase tran- 
sition. 

(ii) In several heavy-fermion metal a 11 ! 12 , concurrent di- 
vergence of the Sommerfeld ratio 7(T) = C(T)/T 
and the magnetic susceptibility xCO * s observed 
at a point that is separated by an intervening NFL 
phase from termination points of any second-order 
phase transitions. 

(iii) In many instances of well-pronounced NFL behav- 
ior, the order parameters required to specify asso- 
ciated second-order phase transitions are still elu- 
sive, casting further doubt on the fluctuation sce- 
narios. 

(iv) In external magnetic fields, thermodynamic proper- 
ties demonstrate scaling behavior governed specifi- 



cally by the ratio fifH/T where fj,f is the magnetic 
moment of constituent fermions. 

These NFL phenomena can be understood when one 
recognizes that standard FL theory possesses its own 
quantum critical point, in the vicinity of which it fails. 
At this point, the necessary stability condition (NSC) for 
the T = Landau state is violated^ 3 - - — as opposed to 
violation of some PSC at a conventional QCP. 

The NSC states that an arbitrary admissible variation 
6n(p) from the FL quasiparticle momentum distribution 
n F(p) — 6(pf — p)i while conserving particle number, 
must produce a positive change of the ground-state en- 
ergy E , i.e., 



SEq = / e{p;np(p))8n{p)dv > 0. 



(1) 



Here, e(p; up) denotes the spectrum of single-particle ex- 
citations measured from the chemical potential [a(T = 0) 
and evaluated for the initial Landau state specified by the 
quasiparticle occupancy np(p). The reduction in energy 
due to breakdown of the NSC, which involves contribu- 
tions linear in <5n, is clearly larger than that due to vio- 
lation of any PSC, which involves bilinear combinations 
of Sn. We must conclude that any associated fluctuation 
scenario is irrelevant to the different type of QCP asso- 
ciated with violation of the NSC, which we shall call a 
Fermi-liquid QCP. 

Violation of the NSC (H} is unambiguously linked to a 
change of the number of roots of equation 



e{p,n F ) = 0. 



(2) 



In standard Fermi liquids, this equation has a single root 
at the Fermi momentum p p , and in that case the signs of 
e(p) and Sn(p) coincide, ensuring satisfaction of the NSC 
([1]). However, consideration of the full Lifshitz phase di- 
agram anticipates the emergence of additional roots of 
Eq. ©. For example, such roots appear at a critical 
density p<> where the function e(p, p ) attains either a 
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The difference in critical indexes is associated with the 



FIG. 1: Three scenarios of emergent bifurcation in Eq. p|: 
Pb < Pf (top panel), pt > Pf (middle panel), pb = Pf (bottom 
panel) . 



maximum, with a bifurcation point pb < Pf, or a mini- 
mum, with pi, > pf, so that e(p — » Pb,Po) <x (p — Pb) 2 , 
(see the upper two panels of Fig. I). Thus, vanishing of 
e(j>b, Po) is always accompanied by vanishing of the group 
velocity v(pb,p<>) — (de(p, p (> )/dp) Pb . Beyond the critical 
density p , the NSC fails to hold, since e(p,riF',p) and 
dn(p) have opposite signs close to pb- 

As indicated in the lower panel of Fig. 1, the condi- 
tion (jTJ) is also violated at a critical density poo where 
the effective mass M*(p) diverges. In this case, standard 
manipulations based on the Landau relation connecting 
the single-particle spectrum and the quasiparticle mo- 
mentum distribution, (see Eq. (J9j) below) yield 



MP) = M _ -i _ 1 po fn \ 



(3) 



where v F 



Pf/M and Ff(fi) = fi(p F ,p F ; p)pfM/it 2 
is the dimensionless first harmonic of the Landau inter- 
action function, normalized with the density of states 
No = pfM/tt 2 of the ideal Fermi gas. Evidently, F°(p) 
is a smooth function of the density p, and F®(p) = 3 
at p = poo. Then beyond the critical point, one has 
-Fi(p) > 3, and the Fermi velocity vf{p) becomes nega- 
tive. This behavior conflicts with the fluctuation scenario 
for the QCP, in which such a sign change is impossible. 

To summarize, we infer that at any point where the 
NSC is violated, the density of states, given by 



AT(T) = I / „( ( ,(1 



(4) 



diverges at T — ¥ due to vanishing of the group velocity 
de(p)/dp. One has^iiZ 

N(T — > 0, poo) oc T -2 / 3 



N{T^0,p o ) ocT- 1 / 2 . (5) 



fact that dp/de 



oc e 



-2/3 



at the critical density p a 



whereas dp/de oc |e| -1 / 2 at the critical density p 

Significantly, the Sommerfeld- Wilson ratio Rsw — 
x{T)/j(T) cannot diverge at these points. Indeed, the 
density of states N(T) cancels out in the ratio Rsw, 
while the Stoner factor entering x(T) maintains a finite 
value, since, as we have seen, the PSC and NSC cannot 
fail at the same point. This conclusion is in agreement 
with experimental data £ 18 ' 19 on dense films of liquid 3 He, 
the two-dimensional electron gas of MOSFETs, and the 
majority of heavy-fermion metals. 

Since no symmetry is violated at a Fermi liquid QCP, 
and hence no hidden order parameters are involved, the 
transition ensuing from the violation of the NSC (p} 
is topological in character! 20 ' 21 Beyond the bifurcation 
point, Eq. ([2]) usually has two additional roots m and 
pi situated near each other (however, cf. Refs. l22H24[ ). It 
is for variations 6n(p) involving momenta p\ < p < p2, at 
which Sn(p) and e(p) have opposite signs, that the NSC 
(fTJ) breaks down. 

The analysis of topological rearrangements triggered 
by the interaction between quasiparticles began twenty 
years agoj 2 ^ with important subsequent developments re- 
ported in Refs. EUa-Hl 

In this article, we address the 
Fermi-liquid QCP in homogeneous matter and focus on 
the case where the new roots p\ and P2 emerge near the 
Fermi momentum pp. The physics of this phenomenon 
is captured if we keep the three first terms, 

e(x) = p F x [v F + — x + — x J, 

v(x) = vf + vix+^x 2 , (6) 

in the Taylor expansions of the spectrum e[x) and its 
group velocity v(x), where x — (p — pf)/pf- To some 
extent, this approach is reminiscent of that employed by 
Landau in his theory of second-order phase transitions. 
In an ideal Fermi gas, v F = v% = v F = (2Me%) 1 ^ 2 . The 
case v\ = 0, V2 > was considered in Ref. [lT]. Here we 
assume that V\ > 0, t>2 > 0, and V1/V2 -C 1, the situation 
addressed in the numerical calculations of Ref. 

To find the bifurcation momentum pb — pf(1 + ^b) one 
must solve the set of equations e(p) = and v{p) = 0, 



V F + -Z-Xb + —X b 

2 o 

, V2 2 
V F + ViX b + —X b 



This system has the solution Xb 
critical condition 

8v 2 v F (p) 



0. 
0. 



(7) 



3« 2 



-3t>i/2t>2 provided the 



(8) 



is met. Thus in the case V\ ^ 0, the critical Fermi ve- 
locity Vf is still positive, and therefore the Landau state 
becomes unstable before the system reaches the point at 



3 



which the effective mass diverges — as was first discovered 
and discussed in Refs. [25|. 

The prerequisite xt <C 1 for applicability of the expan- 
sion ([5]) is satisfied provided v\jvi <C 1, implying that 
the critical Fermi velocity is small: vp — 3v 2 /8v2 -C v F . 
Given this situation, upon accounting for the depen- 
dence of vf on the temperature T and control param- 
eters such as the external magnetic field H that do not 
change the form of Eq. ((8|), one can establish a critical 
line T = T (H) separating phases with different topolog- 
ical structure. 

Evaluation of relevant T— and if— dependent correc- 
tions to the Fermi velocity vf is based on the Landau 
equatio n 33 ' 34 for the single-particle spectrum e(p), which 
in 3D has the form 



de(p) 
dp 



P + l_ 
M 3 



dn( Pl ) 

hip, Pi) — s dv 1 , 

dpi 



(9) 



with dv — p 2 dp I it 1 . This relation provides a nonlin- 
ear integral equation for self-consistent determination of 
e(p,T,H) and the momentum distribution 



n{p,T,H) 



l + e 



■(p,T,H)/T 



(10) 



with the Landau interaction function f(p, pi) (hence its 
first harmonic /i) treated as phenomenological input. 

Our goal is to evaluate the T— and H— dependence 
of the key quantity vp(p,T, H). In the simplest case 
H = 0, the overwhelming T— dependent contributions to 
vf come from integration over the vicinity of the bifur- 
cation momentum pb- Evaluation is performed along the 
same lines as in Ref. Il7l . i.e., by expanding the interac- 
tion function in a Taylor series, although here we have to 
retain a correction to the FL formula Q linear in p — pb- 
As a result, we arrive at 

v F (T -> Tc,) - v F (p) oc J (s-p b ) dn ^ s T) ds, (11) 

where vf{p) is given by Eq. ([3]). 

The integral It on the right side of Eq. (fTTj) is evalu- 
ated with the aid of relations e(p — > pb) oc (p — Pb) 2 and 
de(p — > pb)/dp oc \J e(p) stemming from Eq. ([6]). Upon 
standard changes of integration variables p — > e — > Tz, 
we find 

It oc T 1/2 J z 1/2 n{z) (1 - n{z)) dz oc ^T/e° F (12) 

at T — >• T . Together with Eq. (0), this result leads to a 
tiny value of the critical temperature 



oc e. 



V2Vp 



D 



(13) 
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FIG. 2: Occupation numbers n(p) (top panel), single-particle 
spectrum e(p) in units of 10~ 3 €% (middle panel), and ratio 
e(p)/T (bottom panel) evaluated for the model (|24|) with k = 
0.07 and g s = 0.45, at four color-coded temperatures (in units 



of e%) below T* = 3 x 10~ 



3 



is a criticality parameter. Temperature T<> vanishes at 
the point where D = 0. 

Consider now the imposition of a magnetic field H on 
the system. The impact of the field becomes well pro- 
nounced when PfH > T, and there emerge two subsys- 
tems having spin projections ±1/2, implying in turn a 
decomposition N(T) = N + (T) + N-(T) of the density of 
states. The corresponding formulas are cumbersome and 
will be analyzed elsewhere. Here we focus on the case 
T = and estimate an upper tuning magnetic field H 
such that a bifurcation emerges in the spectrum e+(p), 
while the down-spin spectrum e_(p) = e+(p) — 2/j.fH 
admits merely the conventional root p~p. As before, a 
leading correction Ih to vf comes from integration over 
the vicinity of the momentum pb, with the subsystem 
whose spectrum goes to e+(p) oc (p — p^) 2 — Tz + iifH 
at H — > _ff making the dominant contribution, to yield 



where 



and 



H oc / (Tz + \ijHv) 112 n{z) (1 - n(z)) dz oc ^fj,fH<>/e F 

(15) 



D = 1 — 



3v 2 



1 - 3-Ff CP) ) >0 



(14) 



HfHo oc e° F 



D 



(16) 
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Comparing Eqs. (IT3l) and (TT6]) we see that /ifH ~ T . 
This result is inherent to a scenario in which single- 
particle degrees of freedom play the dominant role and 
is consistent with available experimental data on heavy- 
fermion metals i 19 i 35 

Let us now briefly analyze the situation at T = H = 
on the ordered side of the topological rearrangement as- 
suming, as before, the criticality parameter D to be pos- 
itive. In the case v\ > 0, addressed first in Refs. [H and 
later in Ref. flH the bifurcation momentum p^ resides 
inside the Fermi volume. The rearranged T = quasi- 
particle momentum distribution nip) is given by nip) = 1 
for p < pi and pi <p < pf, and zero otherwise, with pp 
shifted outward to conserve quasiparticle number. Thus, 
the Fermi surface gains an additional hole pocket. In 
the 1960's, such a small hole pocket was called a Lifshitz 
bubble (LB) in Landau-school folklore. In this case, two 
additional roots of Eq. ^ appear, with 



Zl,2 




(17) 



We note that vpip,Pi,P2) differs from the parameter 
vf(p) introduced previously, since it is evaluated for the 
phase in which the Fermi surface has three sheets. Ac- 
counting for the displacement of pf due to emergence of 
the LB, one obtains vp{p,Pi,P2) — v f{p) oc (pi — P2) 2 , 
leading to 



P2 - Pi oc —VD. 

V-2 



As a result, we find 



v i 1 — 
VLB oc — V D < v F 
v 2 



(18) 



(19) 



for the LB Fermi velocity vlb — v(xi) from the second of 
Eqs. (O, thereby demonstrating that the LB contribution 
to the density of states N(0) prevails. 

At temperatures beyond T > T<>, the LB contribution 
to thermodynamic properties disappears. Were this to 
occur instantaneously, the specific heat C(T) would un- 
dergo a jump, as if one were dealing with a second-order 
phase transition. As a matter of fact, the rearrangement 
occurs rapidly but not momentarily. Thus one deals with 
a topological crossover (TC), and Eqs. ([7]) serve to estab- 
lish a TC line T (H) that resembles a line T N {H) of 
second-order phase transitions. 

The TC width is found from the condition T< < T < 
T> , with the boundaries T< and T> being determined by 
the relations 



Since v(p,T) 



e(Pfc,7>) = -T 



(20) 



Ve(p) \JT l<P F in the LB region near T , 
the similar formula e(p,T) = (e<>(p) — eo(pF))VT is ob- 
tained for the spectrum e(p, T) after a simple momen- 
tum integration. Straightforward manipulations employ- 
ing the definition e(pb,T ) = then lead to 

T> - T T - T< 



To 
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FIG. 3: Same as in Fig. [2] but at T > T* 



Accordingly, the reduced temperature width of the crit- 
ical region turns out to be small, implying that the TC 
does indeed imitate a second-order phase transition. 

A conventional FL regime having T-independent quan- 
tities x(0) oc 7(0) cx N(0) (x l/v LB (0) ^ 1/[v f (0)VD] 
is seen to persist until T reaches T< < T , where the 
LB occupation numbers begin to experience substan- 
tial change as the temperature continues to increase. 
Both the density of states N(T) and the spin suscep- 
tibility x(T) attain maximum values at T — T>, where 
X(T>) oc l/(v(xb, To)) cx l/vT©. At higher temperatures, 
the LB contribution to x(T) begins to fall, finally dying 
out and leaving x{T) oc 1/^(0). Analogous results are 
found for the Sommcrfcld ratio, given by 



7(T) = 



t{p) dnjp) 
T dT 



dv, 



(22) 



(21) 



except that 7 reached its maximum at a different tem- 
perature, due to the marked dependence of the spectrum 
e(p, T) on T. The foregoing analysis therefore leads to the 
conclusion that in the QCP region, both the magnetic 
susceptibility x(T) and the Sommerfeld ratio 7(T) ex- 
hibit asymmetric peaks, located at different temperatures 
~ T<>. Such behavior, observed in many heavy-fermion 
metals situated in a QCP region; 19 ' 35 i 36 remains unex- 
plained within any conventional scenario for the QCP. 

Further temperature evolution of the spectrum e(p, T) 
is associated with another essential rearrangement^ of 
the momentum distribution n(p,T) that occurs in the 
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FIG. 4: Same as in Fig. [3] but for the interaction function 
(|25[l with the parameters a = 10 and g y = 70 at three line- 
type-coded temperatures. 



region of a critical temperature The distribution 

T) becomes a smooth function of momentum 



n(p, T) ~ n* (p) , Pi <p < p f , 



(23) 



in an interval adjacent to the Fermi surface and is other- 
wise unity for p < pi and zero for p > pf. In this domain, 
n(p, T) is nearly independent of T, while the dispersion 
of the single-particle spectrum e(p, T) becomes propor- 
tional to T so as to satisfy Eq. (|TU| . 

Both these features are inherent to the phenomenon of 
fermion condensation, a topological phase transition dis- 
covered twenty years agO ) 13 ' 22 ~ 24 i 37 ' 38 in which a fiat band 
pinned to the Fermi surface (the so-called fermion con- 
densate (FC)) is formed. This phenomenon, alternatively 
viewed as a swelling of the Fermi surface, was recently re- 
discovered by Lee 3 ^ while investigating the finite-charge- 
density sector of conformal field theory (CFT) within 
the AdS/CFT gravity /gauge duality. The phenomenon 
of fermion condensation (fiat band) may also arise in 
topological media for purely topological reasons, (see 
RefsliqlilT). 

Unfortunately, important details of this rearrangement 
cannot be established analytically. To clarify the rela- 
tionship between properties of the phase having the sin- 
gle LB at T — and those of a system possessing a 
FC at T — 0, we must resort to numerical treatment of 
Eq. ©. Figs. [3] and |3] present results from numerical 
calculations^ of the spectra e(p) and momentum distri- 
butions n(p) for a 3D model system based on the inter- 




FIG. 5: Spin susceptibility x(T) (top panel) and Sommer- 
feld ratio 7(T) = C(T)/T (bottom panel), in Fermi gas units 
PfM/tv 2 and pfM/3, respectively. Red curves show results 
for the model (1241) and blue curves, for the model (1251). 



action function 



f(q) = 9. 



n 2 p F 



M 



(24) 



with dimensionless parameters g s — 0.45 and /3 = 0.07, 
values for which the zero-T phase possesses a single LB. 



In this model one has T a 



5 x 10- 5 e° P and T* 



10" 



The results are to be compared with those in 



Fig. [4] obtained for the model interaction function 15 



/(?) = 9y 



-aq/pF 



M 



(25) 



for which a flat portion in the spectrum e(j>) is already 
present at T = 0. In the interval T ~ T < T», the spec- 
tra e(p) of the two systems are quite dissimilar. How- 
ever, when T reaches values around T*, a flat portion 
of e(p) develops for the interaction model (|2~4"|) as well, 
the density associated with the flat segment being half 
the FC density obtained for the model (|251) . On the 
other hand, outside the range [pi,pf], the momentum 
distribution n(p) calculated for model (|24|) shows more 
pronounced tails than in the case of model (|2"5j) . 

The impact of these differences on the magnetic sus- 
ceptibilities x an d Sommerfeld ratios ~/(T) of the two 
model systems is seen in Fig. [5] As is known^M 3 . the 
contribution \* of the FC region to x obeys the Curie 
law: x*(T) oc C*/T, with an effective Curie constant 



C* oc / n»(p)(l — n*{p))dv oc p* 



(26) 



Fig. [5] shows that this Curie- like term does prevail in the 
susceptibilities calculated for both models. Such NFL 
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FIG. 6: DC susceptibility (in Fermi gas units) as a function of 
temperature (in units of s%) for different external-field magni- 
tudes. Colors correspond to different values of jj,fH/e%: 10 -3 
(red curve), 3 x I0~ 3 (blue curve), 5 x 10~ 3 (brown curve), 
9 x 10 -3 (green curve). Top panel: model (|24[) . bottom panel: 
model ((231) . 



behavior is in agreement with experimental data on dense 
films of liquid 3 He obtained in the QCP regioni 7 ' 8 ' 10 . 

The Sommerfeld ratios Rsw = x{T)h{T) evaluated 
for the two models differ drastically. Indeed, in the LB 
case (i.e., for model (|24p). the values of \(T — > 0) and 
7(T — > 0) are large compared with those of the corre- 
sponding ideal Fermi gas, with Rsw remaining of order of 
unity. On the other hand, in the model (|25p that already 
hosts a FC at T = 0, x(T) diverges at T -> as T" 1 , 
while the FC contribution to j(T) evidently vanishes, im- 
plying a huge enhancement of Rsw- This scenario is in 
agreement with data^ on the specific heat C(T) of the 
P-type heavy- fermion metal YbIr2Si2, for which the FL 
term C(T) cx T exists only at extremely low tempera- 
tures below To, whose value is presumably less than 1 K. 
The corresponding value of j(T — > 0) is so enhanced that 
already at T ~ 0.7 K, the entropy value has become sur- 
prisingly large: S/N ~ (ln2)/2. At T > T<> a collapse of 
C(T) occurs, and the Landau term linear in T completely 
disappears. Thus, the P-type of the compound YbIr2Si2 



presents the first example of a new class of metals, in 
which flattening of the single-particle spectrum results 
in an ordinary NFL shape of C(T) at extremely low T. 

Next, recall that in conventional Fermi liquids, the 
value of the susceptibility x(T, H) is proportional to the 
density of states at finite field H, but is almost indepen- 
dent of H . By contrast, as witnessed in Fig. [51 an exter- 
nal magnetic field suppresses the NFL contribution to x 
in Fermi systems whose spectra e(p) exhibit a flat region 
(and to the same extent independent of which model is 
chosen). The larger the magnitude of the dimensionless 
parameter HfH/T, the more pronounced is the suppres- 
sion. These results, evaluated within the scheme elab- 
orated in Refs. Il7ll3lll44l . elucidate the NFL behavior of 
dense liquid- 3 He films reported in Refs. FtT-TToI. Taking for 
H a typical value of I T and for \i f the magnitude of the 
3 He atom's magnetic moment, the inequality /i/P > T 
is met at T < 0.5 mK. Since temperatures below 0.2 mK 
are currently attainable, we suggest that it is experimen- 
tally feasible to verify or refute the predicted existence of 
a domain of the (T, H) phase diagram of 2D liquid 3 He 
sensitive to the magnitude of H. 

We have investigated topological transitions arising in 
strongly correlated Fermi systems beyond a quantum 
critical point of the Fermi-liquid type, at which the den- 
sity of states diverges while the Sommerfeld- Wilson ratio 
remains finite. We have attributed these transitions to 
violation of the necessary condition for stability of the 
Landau state. We have shown that in the QCP den- 
sity region, the relevant phase diagram features different 
topological crossovers, occurring between states of the 
same symmetry but with different numbers of sheets of 
the Fermi surface. Importantly, the QCP scenario based 
on topological crossovers does not entail any order pa- 
rameters; hence it is free from the persistent ambiguity of 
conventional fluctuation scenarios associated with hidden 
order parameters. Our analysis predicts the existence of 
a domain of the (T, H) phase diagram of 2D liquid 3 He 
that is sensitive to the magnitude of the magnetic field. 
This prediction is subject to experimental test. 
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